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Power spectral density measurements of any sampled signal are typically restricted by both acquisition 
rate and frequency response limitations of instruments, which can be particularly prohibitive for video-based 
measurements. We have developed a new method called Intensity Modulation Spectral Analysis (IMSA) that 
circumvents these limitations, dramatically extending the effective detection bandwidth. We demonstrate this by 
video-tracking an optically-trapped microsphere while oscillating an LED illumination source. This approach 
allows us to quantify fluctuations of the microsphere at frequencies over 10 times higher than the Nyquist 
frequency, mimicking a significantly higher frame rate. 
© 2009 Optical Society of America 
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Measuring the power spectral density (PSD) is a useful 
way to characterize fluctuations and noise for a diverse 
range of physical processes. Optical measurements of the 
PSD have been used to study single molecule dynam- 
ics [1], bacterial chemotaxis and motion [2,3], quantum 
dot blinking [4] and micror neology [5,6], and to calibrate 
optical traps [7]. Unfortunately, limited acquisition rate 
and detector frequency response restrict PSD measure- 
ments in frequency space, which is especially detrimen- 
tal to video applications. Notably, the highest frequency 
that can be sampled directly, with few exceptions [8-11], 
is half the acquisition rate, or Nyquist frequency. 

We have developed a method called Intensity Modu- 
lation Spectral Analysis (IMSA) , which overcomes these 
limitations in a simple and economical way. By simply 
oscillating the light intensity of an optical signal prior 
to detection, the PSD at the oscillation frequency can 
be determined from the measured variance, even if that 
frequency is above the acquisition rate. This is similar in 
spirit to signal processing methods that can spectrally 
shift a signal (e.g. heterodyne detection, lock-in tech- 
niques), or extract high-frequency information folded 
down via aliasing (e.g. undersampling [11,12]). Practi- 
cally, IMSA can dramatically extend the frequency range 
of an existing measurement device, allowing, for exam- 
ple, an inexpensive camera to be used in place of a signif- 
icantly more expensive one. Here we present the frame- 
work for IMSA and an experimental demonstration using 
an optically-trapped microsphere, in which modulating 
the brightness of an LED allows PSD measurement well 
beyond the Nyquist frequency and camera frame rate. 

IMSA: Concept and foundations Physical acquisi- 
tion systems do not make instantaneous measurements 
when sampling a signal, but rather collect data over finite 
integration times. Consider a stationary random trajec- 
tory X(t). The measured trajectory X m (t) can be ex- 



pressed as a convolution of the true trajectory and an 
impulse response H(t): 

X m (t) = (X * H){t) ee / X(t')H(t - t')dt' (1) 



Then the measured power spectrum P m (oj) differs from 
P(to), the true power spectrum of X, according to the 
relation P m (a>) = P(lu)\H(uj)\ 2 , and the total measured 
variance is the integral of P m (uj) over all frequencies: 
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where the tilde designates the Fourier transform, X(u>) — 
J X(t) cxp(— iujt)dt . Integrals are taken from — oo to 
+oo unless otherwise specified. 

In the simplest case, the measured quantity X m is the 
unweighted time average of the true value X over the 
integration time W, i.e. the impulse response is a rect- 
angular function: 
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elsewhere 
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Correspondingly, the measured power spectrum is the 
original power spectrum multiplied by: 



sin(u>W/2)Y 
ujW/2 ) 
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This simple case of a rectangular impulse response is a 
good model for video- imaging acquisition systems, where 
W is the exposure time. This averaging leads to the 
common problem of video image blur, which not only 
adds errors in the position of tracked objects, but also 
causes systematic biases when quantifying fluctuations 
[10,13-15]. As we have previously demonstrated [10], by 
measuring the variance for different exposure times W 
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and fitting to equation 2, the power spectrum can be 
characterized above the acquisition rate of the detection 
system. However, this approach requires that the func- 
tional form of the power spectrum is known a priori. 

Interestingly, by oscillating the intensity of a source 
signal (e.g. light in the case of video imaging) and 
measuring the variance of the resulting signal, the power 
spectrum can be reconstructed without prior knowledge. 
This is the fundamental idea behind IMSA. 

During detection, the finite integration time W causes 
the true power spectrum to be multiplied by the low- 
pass filter |ifo(w)| 2 (equation 4). As W becomes longer, 
|-£/o(w)| 2 approaches an unnormalized delta function, i.e. 
\H (oj)\ 2 — » aS(uj). Multiplying the rectangular impulse 
response by a complex exponential shifts this delta func- 
tion, i.e. if H(t) = exp(-iu't)H (t), then \H{uj)\ 2 -> 
a6(u> — uj'). From equation 2 we see that the variance 
approaches aP(uj')/2w. Thus, as demonstrated in this 
simple example, the power spectrum can be sampled by 
shifting the filter to any frequency of interest u' and 
measuring the variance. 

Practically, we use the real-valued impulse response 



H(t) = L(t)H (t)/N 



(5) 



where L(t) = sin(a/i + <f>) + B and A is a normalization 
factor (sec Fig. 1). We set N = f L(t)H (t)dt so that the 
convolution represents a time-averaged signal weighted 
by the source intensity modulation L(t), and let B > 0: 
the parameters (W, B) should be chosen such that N ^ 
0. Note that N = B whenever <f> (the phase-shift between 
modulation and sampling) is an integer multiple of 27r, 
or whenever there are an integer number of oscillations 
within each exposure window. 

By substituting H(uj) (the Fourier transform of equa- 
tion 5) into equation 2, and noting that P{—lu) = P{oj) 
for a real- valued wide-sense stationary process, we obtain 
for the measured variance: 
1 
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This key equation can be used to determine the power 
spectrum from the measured variance, since term tl ap- 
proaches P(u')/2WN 2 as W gets longer (see endnote 
[16] and Fig. 1), while terms tS and H become small and 
term t2 can be measured directly. Additionally, proper 
phase selection can completely eliminate t3 and t4- 

IMSA Usage To determine the power spectral den- 
sity of a signal at the frequency of interest a/, the sig- 
nal should be convolved with L(t)Ho(t)/N (equation 
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Fig. 1. Graphical representation of IMSA showing (a) 
unnormalized impulse responses H(t) = L(t)Ho(t) for 
= and B = (solid line, see equation 5) and H(t) = 
H (t) (dotted line, see equation 3), and (b) the squared 
magnitude of their Fourier transforms, which filter P{uj) 
in terms tl and t2 of equation 6. 



5) by modulating its intensity, and then its variance 
should be measured. Writing the measured variance as 
var[X m ](a/), the PSD is given by the following formulas: 
When B = 0, 



P{uj') = 2WN 2 ™[X m ]{uo') 

provided <fi = 7r/4 + m7r/2, where m is an integer. 
When B^O, 
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P(u') = 2WB 2 (var[A m ] (J) - var[A m ] (0)) (8) 

provided that the exposure time W is chosen such that 
W = n2n:/u}' where n is a natural number, and cj> is 
chosen according to (see endnote [17]): 



(-l^shr 1 (sjB 2 + \j2-B 



(9) 



On a practical note, terms tS and H are small when 
W ^> l/u>' (i.e. there are many oscillation cycles within 
each exposure window). If both cross-terms t3 and t4 
are negligible, the selection of 4> is irrelevant. In fact, os- 
cillations may not even have to be synchronized to the 
acquisition device to make a good IMSA measurement — 
simply multiplying the input signal by L(t) prior to de- 
tection can yield acceptable results. 

The error in P(oj') is governed primarily by the error 
in the variance (e.g. for N samples the relative statisti- 
cal standard error is y/(2/N), barring instrumental er- 
ror [10]). The resolution in uj' is given by the width of 
|i/o(o;)| 2 (see Fig. 1), which is approximately ir/W in 
each direction (77% of the area under the curve). 

Experimental Demonstration To demonstrate 
IMSA experimentally, we measured the power spec- 
trum of an optically trapped polystyrene microsphere 
(2.5 um, Corpuscular) using a machine vision camera 
(GE680, Prosilica) and an LED (MRMLED, Thorlabs) 
with intensity modulation capable of up to 300kHz 
(custom current source, Rowland Institute electronics 
lab) . Details of the particle tracking have been discussed 
previously for a functionally identical setup, where it 
was demonstrated that the system is well-described by 
equation 3 [10]. The light intensity was sinusoidally 
modulated (as in equation 5) with B = 0.5 and <fi 
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Frequency [Hz] 

Fig. 2. Measurement of the power spectrum using IMS A 
for two trapped beads with different spring constants 
(data points). The expected power spectra are superim- 
posed (solid lines), showing good agreement well beyond 
the Nyquist frequency (vertical dashed line). Error bars 
represent statistical error. 



determined from equation 9, and the integration time 
of the camera was W = 2.5 ms. Modulated light at 5 
different frequencies lu' was interspersed with DC light 
on a frame by frame basis. Each of the 10 variances 
was calculated and each PSD data point was calculated 
according to equation 8, where the mean of the two 
neighboring DC frames was used as var[X TO ](0) for 
each of the 5 measurements. Data was collected for two 
beads at different laser powers. The spring constant and 
friction factor for each bead was measured using the 
blur-corrected power spectrum method [10]. The IMS A 
measured values, which have no free fitting parameters, 
arc in excellent agreement with the expected power 
spectra (Fig. 2). 

Concluding Discussion As we have shown, IMSA 
provides a unique and practical way for measuring the 
PSD of a signal, which overcomes acquisition rate and 
frequency response limitations of instruments. While es- 
pecially beneficial to video imaging applications (due to 
the almost negligible cost of implementation compared 
with the expense of fast video cameras), IMSA is a gen- 
eral method applicable to signals in which the time- 
averaged weighting can be controlled before sampling. 

Owing to its potential to benefit a broad range of fields 
and its flexibility of implementation (see endnote [18]), 
we expect that IMSA will become a common labora- 
tory method for measuring power spectra. Microrheol- 
ogy measurements [6,19] can take immediate advantage 
of the up to ~ 10000 fold increase in frequency range 
over standard video imaging, without losing the abil- 
ity to track multiple targets. Widely used fluorescence- 
imaging techniques [20] stand to benefit as well, with 
IMSA enabling measurements at frequencies which are 



currently impossible by any other method owing to light 
limitations. 

The dramatic increase in frequency range enabled by 
IMSA can be used to push the envelope of high frequency 
measurement and to realize significant cost savings in 
instrumentation. Using an inexpensive LED illuminator, 
we transformed a standard video camera into a spectrum 
analyzer with a frequency range of up to ~ 300 kHz. 
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